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ABSTRACT 

A 5'x5' region encompassing NGC 63341 and I(N) has been mapped at a wavelength 
of 3 mm (from 83.5 to 115.5 GHz) with the Mopra telescope at an angular resolution 
between 33" and 36". This investigation has made use of the recently installed 3 mm 
MMIC receiver and the Mopra Spectrometer (MOPS) with broadband capabilities 
permitting total coverage of the entire frequency range with just five different ob- 
servations. In total, the spatial distribution of nineteen different molecules, ions and 
radicals, along with additional selected isotoplogues have been studied. Whilst most 
species trace the sites of star formation, CH3CN appears to be most closely associated 
with NGC 63341 and I(N). Both CN and C2H appear to be widespread, tracing gas 
that is not associated with active star formation. Both N2H + and HC3N closely resem- 
ble dust continuum emission, showing they are reliable tracers of dense material, as 
well as the youngest stages of high mass star formation. Hot (E u /k >100K) thermal 
CH3OH emission is preferentially found towards NGC 6334 1, contrasting with I(N), 
where only cold (E u /k <22K) thermal CH3OH emission is found. 

Key words: circumstellar matter - infrared: stars. 



1 INTRODUCTION 

NGC 6334 is a Riant m olecular cloud located at a distance 
of 1.7 kpc l|Neckellll978l ) in the southern Galactic plane. Ly- 
ing along a gas filament of llpc in length, NGC 6334 ex- 
hibits several luminous si tes of massive star fo rmation - as 
seen in the far-infrared ( McBrec n et al.l Il979l . five sources 
denot ed using roman nume rals I to V) and radio contin- 
uum ^Rodriguez et alJll982t six sources denoted with let- 
ters from A to F) with source I dom inating the millime- 
tre to the far infrared (Sandell 2000). In close proximity 
to NGC 6334 1, about two arcminutes to the north, another 
bright source is found in the millimet re- and submillime- 
tre continuum, denoted NG C 6334 1(N) (jCheung et alJll978t 
iGezaril [l982l ; ISandelll lioool ). While this source is believed 
to be a comparably young object - as concluded from the 
lack of Hn regions and mid-IR emission - there is com- 
pelling evi dence however, that star fo rmation is going on 
there fe.g. IMegeath fc Tieftrunklll999h . It is this twin core 
system, NGC 6334 1 and I(N), that offers the rare opportu- 
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nity to study the evolution of high mass stars from the same 
parental cloud and in a relatively small spatial region. 

Single dish molec ular line observations sh ow both cores 
to be chemically rich |Thorwirth et all 120031 ) with source I 
being comparable in line density t o prototypical hot cores 
such a s Orion-KL and S g rB2(N ) (Bachill er fc Cernic haro 
19901; IMcCutcheon et ail l200fj| ; iThorwirth et all 120031 ; 
Schilke et alj|2006l h High spatial resolution observations to- 
wards both cores have been conducted: ATCA investiga- 
tions of NH3 emission up to the (6,6) inversio n transition re- 
veals the pr esence of warm gas in both cores (|Beuther et al.l 
120051. I2007T) and SM A continuum observations at 1.3 mm 
l|Hunter et al.ll2u06h resolve each core into a sample of sub- 
cores of several tens of solar masses each, four for source 
I and seven for I(N), demonstrating the formation of star 
clusters. 

Several line mapping studies using standard molecu- 
lar tracers were carried out in the past to obtain informa- 
tion ab out the large scale spatial distribution of molecu - 
lar gas (|Kraemer fc Jackson! Il999l : IMcCutcheon et al.l 2000) 
additionally leading to the detecti on of several molecular 
outflows associated with both cores JBachiller fc Cernicharol 
Il990l : IMegeath fc Tieftrunk|[l999l : iLeurini et alj|2006l ). 
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In the present investigation we have used the Mopra 
telescope for spectral line mapping of a 5' x 5' region around 
NGC 6334 1 and I(N) throughout the entire 3 mm range. 



2 OBSERVATIONS 

Observations were carried out in 2006 on June 25th, using 
the 22m Mopra telescope near Coonabarabran, New South 
Wales, Australia^- We used the recently commissioned 3 mm 
MMIC receiver in combination with the new Mopra spec- 
trometer, MOPS, in broadband mod^EI resulting in an in- 
stantaneous bandwidth of 8 GHz split over four overlapping 
IFs of 2.2 GHz each. At the time of the observations MOPS 
provided 1024 channels per IF, giving a total of 8192 chan- 
nels. The spectral resolution of the observations was 2 MHz 
corresponding to a velocity resolution of about 6kms -1 per 
channel. For NGC 6334, this velocity resolution is about the 
same, or slightly larger than expected line widths, which 
means while we do not expect to spectrally resolve many 
lines, we will not suffer a significant loss in signal due to 
spectral smearing. Thus, this work will necessarily concen- 
trate on morphologies and integrated intensities of emission, 
rather than any analysis of kinematic structure. 

The pointing centre for the maps was chosen to be 
RA(J2000) = 17 20 54.24, Dec(J2000) = -35 46 11.5, which 
is approximately half way between NGC 6334 1 and I(N). 
We used the on-the-fly mapping routine to map an area of 
5' x 5' around this position. Due to incomplete sampling of 
the edges of the map, a slightly smaller area is considered in 
this work. A scanning rate of 3.5 arcsecs -1 was used. Data 
were averaged every 2 seconds which gives an optimum data 
collection rate with minimal smearing of the output. A spac- 
ing of 10 arcseconds was used between each scanning row. 
Assuming a beam of between 33 and 36", depending on the 
frequency l|Ladd et al.ll2005T ). this observing mode resulted 
in a fully sampled map. A single map was obtained in about 
90 minutes, including about 15% of this time for Tsys cal- 
ibration, reference observations and pointing on a nearby 
SiO maser (AH Sco). Poor weather (clouds) affected some 
of the data, which usually appears as faint horizontal stripes 
in the data (see the following sections for details). 



these transitions, we provide a map of the emission, shown 
in Figure [3] 

Figure [2] shows dust continuum emission in the region 
taken using SCUBA on the JCMT (Goran Sandell, private 
communication). These images were taken on June 26th and 
July 1st, 1998, with measured half-power beam widths of 
14.6" at 850/^m and 8" at 450/im. In this Figure, the plus 
symb ols represent the positions of methanol maser emis- 
sion (Walsh et all 1 19981) . In the north, the methanol maser 
site marks the position of NGC 6334 1(N). In the south, two 
maser sites are seen, spaced about 7" apart. Coincident with 
one of these maser sites, and marked by the diamond is the 
UC Hn region NGC 6334 1. The oth er, nearby maser sit e 
marks the position of a hot core (eg. iBeuther et all 120051 ). 
The circle marks the position, and approximate extent of 
the Hn region NGC 6334 E. Strong sources of submm con- 
tinuum emission appear associated with both NGC 6334 1 
and NGC 6334 1(N) . At the position of NGC 6334 1, the dust 
continuum appears to peak at the UC Hn region position, 
rather than the hot core located to the north-west. The con- 
tinuum emission also appears to be e ssentially unresolved 
here. Higher resolution observations bv lHunter et al.l (I2006T ) 
show the continuum emission does fragment into multiple 
sources. 

The continuum emission associated with NGC 6334 1 
appears to be extended over almost an arcminute, corre- 
sponding to a projected length of nearly 0.5 pc. The contin- 
uum emission also appears to peak about 6" to the south- 
east of the maser site associated with NGC 6334 1(N). The 
greyscale continuum images also show the low level dust 
emission, which exhibits a prominent filamentary structure 
running from the south-western corner to the centre-top of 
the field of view. Other, weaker filamentary structures ap- 
pear as a bridge between NGC 6334 1 and I(N), as well as 
possibly connecting NGC 6334 1 and NGC 6334 E, although 
very little continuum emission is coincident with the Hn re- 
gion NGC 6334 E. The very weak filamentary structure ap- 
pears to extend from NGC 6334 1(N) to the north-eastern 
corner of the field of view. Finally, there is some emission 
seen directly to the north of NGC 6334 1(N), which may be 
an extension of the filament seen to bridge the gap between 
NGC 6334 1 and I(N). 



3 RESULTS 

Figure [T] shows spectra at I and I(N) across the 3 mm band. 
Each spectrum was made by integrating the emission over 
a box approximately equal to the size of the beam. The 
spectrum for NGC 6334 1 was centred on the methanol maser 
site that overlaps with the UC Hn region (17 20 53.35, -35 
47 1.6) and the spectrum for NGC 6334 1(N) was centred on 
the methanol maser site (17 20 54.58 -35 45 8.6). Across the 
band, we detect a total of 52 transitions from 19 species. 
Details of each detection are given in Table [T] For each of 

1 For detailed information on the facility visit the Mopra www 
page at http://www.narrabri.atnf.csiro.au/mopra/ 

2 Detailed information on the available receivers and 
backcnds can be found online in the "Technical Sum- 
mary of the Mopra Radiotelescope" available at 
http://www.narrabri.atnf.csiro.au/mopra/mopragu.pdf 



4 DISCUSSION 

4.1 Morphological analysis 

Since the spectral resolution was relatively coarse, we re- 
frain from an analysis of spectrally resolved data-cubes and 
only show and interpret integrated images of the different 
molecular species and isotopologues. Furthermore, the an- 
gular resolution of the Mopra telescope in the 3 mm band 
is ~ 36" (at 90 GHz) corresponding to a linear resolution of 
~65000AU or ~0.3pc. This spatial resolution is not suffi- 
cient to analyze spatial differences within one or the other 
massive star-forming region. However, it allows us to inves- 
tigate the molecular global properties, similarities and dif- 
ferences between the different massive star-forming regions 
present in our field of view. In the following, we will quali- 
tatively describe the global spatial characteristics of the dif- 
ferent molecular species and its isotopologues. 
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Figure 1. Spectra of NGC6334I (lower) and NGC6334I(N) (upper) shown between 83.5 and 91.5 GHz. The frequencies of identified 
lines in each spectrum are indicated with a solid vertical line. Each spectrum was produced by integrating over a square box with each 
side approximately equal to the beam width (36"). Strong lines in the spectrum occasionally produce a ringing artifact close to the line 
frequency, eg. CH3OH in I(N). This is a result of the Gibbs phenomenon and is most pronounced with narrow lines, such as masers. 



H41q: The only covered hydrogen recombination line 
that we could image and is found as expected toward the 
UCH11 region associated with NGC 6334 1. 

C2H: C2H shows a peculiar spatial morphology since 
its main emission peak is associated with NGC 6334 1(N) and 
appears as a weak hole at the position of NGC 6334 E. Fur- 
thermore, extended emission closely follows the dust contin- 
uum ridge shown in Figure [2] This molecule shows no clear 
peak toward the southern source NGC 63341. The dearth of 
emission coincident with both NGC 6334 E and I suggests 



that this molecule is destroyed close to Hn region environ- 
ments. 

Unsaturated hydrocarbons like C2H are known to be 
strong at early evol utionary stages (eg. lMillar fc Neiadlll985l ; 
iTurner et aLll 19991 ) . During the warm-up of the cores it re- 
acts quickly wi th oxygen to form CO through ion-molecul e 
chemistry (eg. ITurner et af] 1 19991 ; iHerbst fc Leund Il986l ). 
These models are consistent with our finding of the C2H 
peak close to NGC 6334 1(N) and much less emission toward 
the hot core NGC 6334 1 or the Hn region NGC 6334 E. Fur- 
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Figure 1 - continued ... Frequency range is between 91.5 and 99.5 GHz. 



thermore. lBeuther et alj (|2008l ) confirm that during ongoing 
evolution C2H abundances decrease toward star formation 
sites. Thus, C2H is a good tracer of gas that is not directly 
associated with active star formation. 

CN: The radical CN appears to be smoothly dis- 
tributed over a region covering NGC 6334 1, E and I(N). One 
transition is slightly stronger at NGC 6334 1, and the other is 
slightly stronger at NGC 6334 1(N). However, these are weak 
peaks, over the extended emission and so it is not possible 
to tell if the difference in relative intensities is real or due to 
noise variations. CN, like C2H, appears to be a good tracer 
of gas not directly associated with active star formation. 



12 CO, 13 CO, C 18 and C 17 0: Carbon monoxide and 
its isotopologues are easily detected in our field of view. 
However, the spatial distribution is surprising because it 
does not trace clearly any of the two main massive star- 
forming regions, NGC 6334 1 or I(N). The main emission 
peak is located approximately 50" north-west of NGC 6334 1, 
closer to the UCH11 region NGC 6334 E (but not coincident). 
We note that this peak does not appear pronounced in any 
other spectral line imaged in this work, nor does it feature 
prominently in the dust continuum map. Since only CO (1- 
0) transitions are prominent at this position (17 20 48.45, 
-35 46 30) and these transitions have a low effective critical 
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Figure 1 - continued ... Frequency range is between 99.5 and 107.5 GHz. The large blank areas shown on the upper plot are because 
the data in these frequency ranges were unfortunately unusable. Note that one line appears at approximately 100.322 GHz. This is a 
reflection of HNC at 90.664 GHz, and not a real line. 



density (~ 10 2 cm -3 ), we believe this line of sight contains 
a high column density of low density gas without any signif- 
icant high density centre that might contain star formation. 
While 12 CO, 13 CO and C 18 show weaker secondary peaks 
toward NGC6334I and I(N), these are barely detectable in 
the rarest i sotopologue C 17 Q. Based on the mm continuum 
maps (e.g., Sandcll 2000, see also Figure[2]) it is obvious that 
the highest gas column densities are toward NGC 6334 1 and 
I(N), therefore, non-detection of these peaks in rare CO iso- 
topologues cannot be explained by column density effects. 
The temperatures toward NGC 6334 1 and I(N) exceed 100 K 



ijBeuther et all l2005f ) whereas the excitation temperatures 
of the upper energy levels of the J = 1 — carbon monox- 
ide transitions are only of the order 5 K. With a typical 
Boltzmann distribution, one would expect, for example, the 
J = 6 — 5 transitions to peak closer to the warm regions 
NGC 6334 1 and I(N). 

An additional interesting feature in the CO maps is the 
more extended north-south ridge. It would be useful to have 
larger maps to investigate to what extend this ridge is con- 
tinuing, but a comparison with the mm continuum map in- 
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Figure 1 — continued ... Frequency range is between 107.5 and 115.5GHz. The frequencies of identified lines in each spectrum are 
indicated with a solid vertical line below the spectrum. 



dicates that this CO ridge has no direct counterpart in the 
dense gas and dust component. 

N2H + : The well-known tracer of star formation at 
young evolutionary stages N2H + exhibits its strongest peak 
toward the northern region NGC6334I(N). However, the 
molecule has additional - although progressively weaker 
- emission peaks associated with the more evolved region 
NGC6334I. N2H + emission is at a local minimum at the 
position of NGC6334E, indicating that this ion tends to 
avoid the Hn region environment. In addition to this, N2H + 
clearly shows extended emission, which is positionally co- 



incident with the main filament shown in dust continuum 
emission (Figure [2} . Of all the spectral line species covered 
in this work, N2H + most closely follows the dust continuum 
emission, showing it to be a good tracer of cold and dense 
material, including in star forming regions. 

HCO+ and H 13 CO+: The HCO+ map clearly shows 
its main emission peaks associated with the dense cores 
NGC6334I and I(N). It shows an additional secondary 
peak associated with the main carbon monoxide peaks fur- 
ther to the east, as well as hints of the larger scale north- 
south ridge discussed for CO and its isotopologues above. 
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Figure 2. Dust continuum images taken with SCUBA on the J CMT, from Goren Sandcll. The large circle represents the approximate size 
and position of the Hn region NGC 6334 E l|Carral et al]|2002t ). The diamond represents the position of the UC Hn region NGC 6334 F 
also k nown as NGC 6334 1. The two plus symbols next to the diamond represent the position of two methanol maser sites i Walsh et al.l 
1998). The single plus symbol in the northern half of each image represents the position of NGC6334I(N), traced by a methanol maser 
site llWalsh et al . 1998). (top) 450/im image. The lowest contour is at 50 Jy beam 1 and contours are spaced at 50 Jy beam 1 . The peak 
flux density in the map is 369 Jybeam -1 . (bottom) 850/^m image. The lowest contour is at 10 Jy beam -1 and contours are spaced at 
10 Jybeam - 1 . The peak flux density in the map is 59 Jybeam - 1 . 



For H CO + we cannot identify the north-south ridge any- 
more but the other features resemble those of the main 
isotopologue. It is interesting to note that the main in- 
tensity peaks for both isotopologues are associated with 
NGC 6334 1(N) and not the more prominent hot-core type 
region NGC 63341. This is rem i niscen t to the HC3N map 
published bv lSollins fc Megeathl (I2004T I. 

CH3OH, thermal and class 1/11 masers: The spec- 
tral setup covers many CH3OH and class I and class 11 maser 
lines as shown in Table [1] (see also lMtiller et al1 [2004). It is 
impossible to definitively claim detection of masers given 
the current observations with a beam larger than 30" , and a 
spectral resolution of a few km/s. However, we can comment 
on the likelihood of any transition showing masing activ- 



ity, based on the relative intensities of lines at the positions 
of I and I(N). It is likely, then, that NGC 6334 1(N) shows 
maser emission in the Class I transitions at 84.521 GHz and 
95.169GHz, and in the Class II transition at 108.894 GHz. 
IVal'tts et al.1 (2000) confirms the masing nature of the Class 
I transition at 95.169 GHz. NGC 6334 1, on the other hand, 
does not show strong emission in any of the masing transi- 
tions, except perha ps th e Class II tran s ition at 107.013 GHz. 
ICragg et al.1 feoOll ) and IVal'tts et af] <| 19991 ) have both ob- 
served this Class II transition towards NGC 6334 1 and con- 
firm it is a maser. 

The thermal methanol transitions appear to divide into 
two types, based on the strength of their emission. Methanol 
transitions showing strong emission at 95.914, the quadru- 
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Figure 3. Integrated intensity imag es and contou rs for detected transitions. The large circle represents the approximate size and position 
of the Hn region NGC6334E dCarral et al.ll2002h . The diamond represents the position of the UC Hn region NGC 6334 F, also know n 
as NGC 6334 1. The two plus symbols next to the diamond represent the position of two methanol maser sites [fWalsh et al.lll998l) . 
The single plus sy mbol in the northern half of each image represents the position of NGC 6334 1(N), traced by a methanol maser site 
fWalsh et aJj|l998ft . Axes are labelled in coordinates of Ra and Dec (J2000). The intensity scale is in units of brightness temperature 
multiplied by velocity, (a) H41o radio recombination line. Contours start at 5a and increase in la steps, where la is 0.9Kkms _1 . 

(b) C2H (N=l— 0, J=3/2-l/2, F=2— 1) — ethynyl radical. Contours start at 5a and increase in la steps, where la is 1.2Kkms _1 . 

(c) C2H (N=l-0, J=3/2— 1/2, F=l— 0) — ethynyl radical. Contours start at 5a and increase in la steps, where la is 1.2Kkms _1 . (d) 
C2H (N=l-0, J=l/2— 1/2, F=l— 1) — ethynyl radical. Contours start at 5a and increase in la steps, where la is 0.9Kkms _1 . (e) CN 
(N=l— 0, J=l/2-l/2, F=2-2) — cyanogen radical. Contours start at 5a and increase in 2a steps, where la is l.lKkms -1 . Note that 
this transition was particularly sensitive to weather. Features in the northern quarter of the image (north of I(N)) are almost certainly a 
result of poor weather and not real, (f) CN (N=l— 0, J=l/2— 1/2, F=2— 1) — cyanogen radical. Contours start at 5a and increase in la 
steps, where la is 1.7Kkms _1 . 
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Figure 3 - continued (g) H 13 CN (1-0) — hydrogen cyanide. Contours start at 5a and increase in la steps, where la is 1.2Kkms 1 . (h) 
HCN (1-0) — hydrogen cyanide. Contours start at 5<r and increase in la steps, where la is 6.0Kkms _1 . (i) HN 13 C (1-0) — hydrogen 
isocyanide. Contours start at 5<r and increase in la steps, where la is 0.6Kkms -1 . (j) HNC (1-0) — hydrogen isocyanidc. Contours 
start at 5a and increase in 2a steps, where la is 2.0Kkms -1 . (k) C ls O (1-0) — carbon monoxide. Contours start at 5a and increase 
in 2a steps, where la is l.lKkms -1 . (1) 13 CO (1-0) — carbon monoxide. Contours start at 5a and increase in 2a steps, where la is 
5-OKkms- 1 . 
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Figure 3 - continued (m) C 17 (1—0) — carbon monoxide. Contours start at 5a and increase in la steps, where la is 0.6Kkms 1 . 
(n) 12 CO (1-0) — carbon monoxide. Contours start at 5a and increase in 3tr steps, where la is 5.0Kkms _1 . (o) N 2 H+ (1-0) - 
diazcnylium. Contours start at 5a and increase in la steps, where la is 5.0Kkms _1 . (p) H 13 CO + (1-0) — oxomethylium. Contours 
start at 5a and increase in 2a steps, where la is OJKkms" 1 . (q) HCO+ (1-0) — oxomethylium. Contours start at 5a and increase in 
la steps, where la is 1.6Kkms — . (r) CH3OH-E (5— 1,5 — 4 q,4) — methanol. C ontours start at 5a and increase in 10c steps, where la 
is 0.6Kkms _1 . This methanol transition is a Class I maser jMuller et alj|2004h . 
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(s) CH 3 OH-A 4 (B 0[8 -7 U ) Class I maser 95169.463 MHz 




(t) CH 3 OH-A + (2 1? -1 1 J 



• '14.309 MHz 




CH 3 OH-A" 



20 m 50 E 30 m 45 : 



^?'o8-i 804 14Hz 




;h 21 n, 05 s 21 m 00 s 20 m 55 s 20 m 50 5 20 m 45 s 



CH 3 0H-E (13 3J1 -12 a9 ) 



lOOe-38.900 MHz 



+ 



(x) CH3OH-E (1L UI -10_ !8 ) Class I maser 104300.414 MHz 



Figure 3 - continued (s) CH3OH-A+ (8q.s— 7 i,7) — methanol. C ontours start at 5<r and increase in ha steps, where la is l.OKkms . 
This methanol transition is a Class I maser llMuller et alj|2004h . (t) CH3OH-A+ (21,2-11,1) — methanol. Contours start at 5<r and 
increase in 2a steps, where la is 0.6Kkms _1 . (u) CH30H-A+ (2o,2 — lo,l) — methanol. Contours start at ba and increase in 5a 
steps, where la is 3.0Kkms _1 . (v) CHsOH-A - (2i,i-li,o) — methanol. Contours start at 5a and increase in 2a steps, where la 
is 0.6Kkms _1 . (w) CH3OH-E (132, 11-1239) — methanol. Contours start at 5a and increase in la steps, where la is 0.3Kkms _1 . 
Strong "emission" is seen on the right hand edge of the image, at a declination of approximately —35° 47', which is due to a bad data 
point, and is not real emission, (x) CH3OH-E (11 — 1.11— 10 — 2.9) — methanol. Contours start at 5a and increase in la steps, where la is 
0.6Kkms — 1 . This methanol transition is a Class I maser (Muller ct al. 20q3). "Emission" is seen on the right hand edge of the image, 
at a declination of approximately — 35°47', which is due to a bad data point, and is not real emission. 
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Figure 3 - continued (y) CH30H-A + (3i.3-4 p,4) — methanol. C ontours start at 5a and increase in 2a steps, where la - is 0.6Kkms . 
This methanol transition is a Class II maser jMuller et al.ll2004T) . (z) CH3OH-E (15_ 2 , 14-151,14) — methanol. Contours start at 5a 
and increase in la steps, where la is 0.6Kkms _1 . Strong "emission" is seen on the right hand edge of the image, at a declination of 
approximately — 35°47', which is due to a bad data point, and is not real emission, (aa) CH3OH-E (Oq.q-1 — 1.1) — methanol. C ontours 
start at 5a and increase in 2a steps, where la is 0.6Kkms _1 . This methanol transition is a Class II maser (Mu ller et al]|2004h . Strong 
"emission" is seen on the right hand edge of the image, at a declination of approximately — 35°46'30", which is due to a bad data 
point, and is not real emission, (ab) CH30H-A+ (72,5— 81, s) — methanol. Contours start at 3<r and increase in la steps, where la is 
0.5Kkms _1 . (ac) CH3CN (5-4) — methyl cyanide. The emission includes all transitions in the K ladder from to 3, inclusive. Contours 
start at 5a and increase in la steps, where la is 0.9Kkms . We attribute the "emission" seen south of — 35°48' to be due to poor 
weather and is thus not real, (ad) CH3CN (6-5) — methyl cyanide. The emission includes all transitions in the K ladder from to 3, 
inclusive. Contours start at 5a and increase in la steps, where la is 0.9Kkms _1 . 
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(ae) HNCO (4 04 -3 03 ) 87925.252 MHz (af) HNCO (5 05 -4 04 ) 109905. 75B MHz 




17 h 21 m 05 5 21 m 00 5 20 m 55 5 20 m 50 5 20 m 45 5 17 h 31 m 05 5 21 m 00 5 20 m 55 5 20 m 50 5 30 m 45 s 



Figure 3 - continued (ae) HNCO (4o,4-3o,3) — isocyanic acid. Contours start at 5<r and increase in la steps, where la is 0.6Kkms _1 . 
Strong "emission" is seen on the right hand edge of the image, at a declination of approximately — 35°47'30", which is due to a bad 
data point, and is not real emission, (af) HNCO (5o,5-4o,4) — isocyanic acid. Contours start at 5a and increase in la steps, where la 
is 0.6Kkms _1 . (ag) 13 CS (2-1) — carbon sulfide. Contours start at 5a and increase in la steps, where la is O.VKkms- 1 . (ah) C 34 S 
(2—1) — carbon sulfide. Contours start at 5a and increase in 2a steps, where la is l.OKkms -1 . Strong "emission" is seen on the right 
hand edge of the image, at a declination of approximately — 35°47'15", which is due to a bad data point, and is not real emission, (ai) 
C 33 S (2-1) — carbon sulfide. Contours start at 5a and increase in la steps, where la is 0.6Kkms _1 . (aj) CS (2—1) — carbon sulfide. 
Contours start at 5a and increase in 5<r steps, where la is 1.8Kkms _1 . 
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17 h 21 m 05 5 21 m 00 5 30 m 55 s 20 m 50 s 20 m 45 s 17 h 21 m 05* 21 m 00 5 20 m 55 5 20 m 50 5 20 m 45 s 



Figure 3 - continued (ak) SiO (2-1) — silicon monoxide. Contours start at 5a and increase in 5a steps, where la is 0.6Kkms _1 . (al) 
NS (5/2-3/2e) — nitrogen sulfide. Contours start at 5a and increase in la steps, where la is 0.6Kkms _1 . (am) H2CS (31,3-21,2) — 
thioformaldchydc. Contours start at 5a and increase in la steps, where la is 1.5Kkms -1 . (an) H2CS (3i,2-2i,i) — thioformaldchydc. 
Contours start at 5a and increase in la steps, where la is 1.3Kkms -1 . (ao) SO (22-I1) — sulfur monoxide. Contours start at 5a and 
increase in la steps, where la is 0.6Kkms _1 . (ap) SO (23-I2) — sulfur monoxide. Contours start at 5a and increase in 5a steps, where 
la is 0.6Kkms — 1 . 
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Figure 3 - continued (aq) SO (32~2i) — sulfur monoxide. Contours start at ha and increase in 1<t steps, where 1<t is 0.6Kkms 1 . 
(ar) HC3N (10-9) — cyanoacetylene. Contours start at 5<r and increase in 2a steps, where la is 1.2 Kleins -1 , (as) HC3N (11-10) - 
cyanoacetylene. Contours start at 5a and increase in 2a steps, where la is l.lKkms -1 . (at) HC3N (12-11) — cyanoacetylene. Contours 
start at 5a and increase in 2a steps, where la is l.lKkms - 1 . (au) OCS (7-6) — carbonyl sulfide. Contours start at 5a and increase 
in la steps, where la is 0.6Kkms _1 . (av) OCS (8-7) — carbonyl sulfide. Contours start at 5a and increase in la steps, where la is 
l.lKkms -1 . Strong "emission" is seen around edges of the bottom half of the image, which is most likely due to poor weather, and is 
not real emission. 
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103463 063 MHz 




(ax) CH 3 OCHO (10 010 -9 09 ) 111171.636 MHz 



17 h 21 m 05 s 21 m 00 s 20 m 55 s 20 m 50 s 20 m 45 : 
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Figure 3 - continued (aw) OCS (9-8) — carbonyl sulfide. Contours start at 5a and increase in la steps, where la is 0.7Kkms _1 . 
"Emission" is seen around the northern edge of the image, which is most likely due to poor weather, and is not real emission, (ax) 
CH3OCHO-A (10o,iq— 9q,9) — methyl formate. Contours start at 5a and increase in la steps, where la is 0.3Kkms _1 . 



plet at 96.7 and at 97.583 GHz all exhibit the bulk of their 
emission from NGC6334I(N). The methanol transitions at 
100.639, 107.160 and 111.290 GHz all show weak emission 
that is only detected towards NGC6334I. 

CH3CN: Methyl cyanide is detected in the rotational 
transitions J — 5 k — 4k and 6 k — 5k- Both integrated 
intensity maps show two comparably strong peaks toward 
NGC6334I and I(N). Methyl cyanide is potentially useful 
for temperature estimates, by comparing intensity r atios of 
components in the if-ladders (|Thorwirth et al. 2003). How- 
ever, the low spectral resolution of these data do not allow 
us to distinguish individual if -ladder elements. Of all the 
species mapped in this work, CH3CN appears to be the one 
that is most concentrated towards the two regions of active 
star formation: NGC6334I and NGC6334I(N). 

12 CS, 13 CS, C 34 S and C 33 S: Carbon sulfide shows 
an interesting spatial variation going from the main isotopo- 
logue CS step by step through the rarer isotopologues 13 CS 
and C 34 S to C 33 S. While the former main isotopologue peaks 
toward the southern hot core region NGC 6334 1, the north- 
ern peak NGC 6334 1(N) is of comparable brightness for 13 CS 
and C 34 S. This is likely due to the main CS isotopologue be- 
ing optically thick, which is then tracing a warmer environ- 
ment around NGC 6334 1. The optically thin tracers 13 CS 
and C 34 S show approximately equal brightnesses in I and 
I(N), suggesting the amount of CS is approximately equal 
in each. However, the rare isotopologue C 33 S is brighter in 
I. We would expect C 33 S to look similar to 13 CS and C 34 S. 
This suprising result needs to be investigated further. 

SO, OCS, H 2 CS and NS: Several sulfur-bearing 
species have been observed, and all of them show two 
peaks associated with NGC 6334 1 and I(N). The general 
trend appears to be that the emission is stronger towards 
NGC 6334 1(N), rather than I. The one exception to this 
rule is OCS (8-7), where very little emission is detected 
towards NGC 6334 1(N), even though NGC 6334 1 shows a 
bright peak. This is an unexpected result, especially in light 



of the fact that the two other OCS transitions clearly show 
NGC 6334 1(N) as a strong source of OCS emission. 

HC 3 N, HNCO, HNC, HN 13 C, HCN, H 13 CN: 
There are several hydrogen-carbon-nitrogen chain molecules 
observed, and all of them show two strong peaks associ- 
ated with NGC 6334 1 and I(N). Except for HCN where the 
two peaks are of approximately the same strength, in all 
other lines the region I(N) is clearly the dominant one. 
For HC3N this has alre ady previously been reported by 
ISollins fc Megeathl (|2004l ). HC3N is unusual in that it shows 
extended emission which appears to follow the main filament 
seen in dust continuum emission (Figure [2]), to the west of 
NGC 6334 E. Furthermore, the HC3N emission appears to 
surround NGC 6334 E, but very little emission is seen in- 
side the ring denoting the extent of this Hn region. Since its 
morphology closely resembles that of the dust continuum 
emission, we find that HC3N is a good tracer of quiescent, 

dense gas, as is the case for N 2 H + . 

SiO: As previously reported by iMegeath fc Tieftrunkl 

(1999), SiO emission is found strongly toward 
NGC 6334 1(N) and is weaker elsewhere in the field of 
view. In NGC 6334 1(N) the SiO emi ssion is clearly associ- 
ated with its molecular outflow(s) (|Megeath fc Tieftrunkl 
1 19991 . Hunter et al. in prep.). Previous observ ations of the 
G333 giant molecular cloud lLo et all <|2007D have found 
SiO emission from a region that shows similar properties 
to NGC 6334 1(N) with an outflow associated with a very 
early stage of high mass star formation. 

CH3OCHO: The dense core tracing molecule methyl 
formate is detected only toward the southern hot core region 
NGC 6334 1. 



4.2 Comparison of NGC 63341 and I(N) 



In i|4.1l we described the morphology of emission for 
molecules, or groups of molecules. Here we synthesise this 
information to present a comparison of the two star forming 
sites NGC 6334 1 and I(N). Based on previous work on this 



Mopra mapping of NGC 6334 1 and I(N) 17 



Table 1. Molecular lines mapped towards t he NGC 6334 1/I(N) regi on in the present study. Transition 
frequencies (MHz) w ere taken from the CDMS jMiiller et a l, 2001, 2005) in most cases and the JPL catalogue 
l|Pickett et al.|[T998h . Also given are the integrated intensities for each transition, as measured at NGC 6334 1 
and I(N) by integrating the emission over one beam, as was done for the spectra presented in Figure [T] 
3cr-upper limits are given for those transitions not detected. 



Molecule Transition Frequency Methanol Integrated Intensity 

(MHz) Maser (K km/s) 

Class NGC 6334 1 NGC 6334 1(N) 



H 








41a 


92034.442 




26 


< 10^ 


C 2 H 




N 


= 1 - 


0, J = 3/2 - 1/2, F = 2 — 1 


87316.898 




66 


60~ 


C 2 H 




N 


= 1 — 


0, J = 3/2 - 1/2, F = 1 — 


87328.585 




38 


32 


C 2 H 




N 


= 1 - 


0, J = 1/2 - 1/2, F = 1 — 1 


87401.989 




35^ 


23 b 


CN 


N 


= 1 


- 0, J 


= 1/2 - 1/2, F = 1/2 - 1/2 


113123.370 




8.9" 


19 


CN 


N 


= 1 


- 0, J 


= 1/2 - 1/2, F = 1/2 — 3/2 


113144.157 




38 


34 


CN 


N 


= 1 


- 0, J 


= 1/2 - 1/2, F = 3/2 — 1/2 


113170.492 




61 


64 


CN 


N 


= 1 


- 0, J 


= 1/2 - 1/2, F = 3/2 — 3/2 


113191.279 




56 


16 


H 13 CN 








J = 1 — 


86339.922 




61 


57 


HCN 








J = 1-0 


88631.602 




350 


300 


HN 13 C 








J = 1-0 


87090.850 




19 


20 


HNC 








J = 1-0 


90663.568 




100 


110 


c 18 o 








J = 1 - 


109782.173 




69 


69 


13 co 








J = 1 — 


110201.354 




360 


350 


c 17 o 








J = 1 - 


112359.284 




19 


21 


12 co 








J = 1 - 


115271.202 




930 


680 


N 2 H+ 








J = 1 - 


93173.392 




210 


440 


H 13 CO+ 








J = 1 - 


86754.288 




16 


24 


HCO+ 








J = 1 — 


89188.525 




170 


150 


CH3OH-E 








JfCa.Kc = 5-1,5 — 4o,4 


84521.169 


I 


80 


180 


CH3OH-A+ 








JK a ,K c = 80,8 — 7l,7 


95169.463 


I 


11 


88 


CH3OH-A+ 








JK a ,K c = 2l, 2 — ll,l 


95914.309 




16 


32 


CH3OH-A+ 








JK a ,K c = 2 0,2 — lo,l 


96741.375 




150 


370 


CH3OH-A- 








Jli a ,K c = — ll,0 


97582.804 




9.5 


24 


CH3OH-E 








JK a ,K c = 132,11 — 123,9 


100638.900 




4.9 


< 0.9' 


CH3OH-.E 








JKo.,K c = H-1,11 — 10-2,9 


104300.414 


I 


5.5 


< 1.8' 


CH3OH-A+ 








Jli a ,K c = 3l,3 — 4o,4 


107013.803 


11 


15 


2.7 


CH3OH-E 








■hi a ,K a = 15-2,14 — 15l,l4 


107159.820 




8.2 


< 1.8' 


CH3OH-E 








JK a ,K c = Oo,0 — 1-1,1 


108893.963 


11 


18 


28 


CH 3 OH-7l+ 








JK a ,K c = 72,5 — 81,8 


111289.550 




1.9 


< 1.2' 


CH3CN 








Jk = 5o - 4o 


91987.091 




41_c 


35 c 


CH 3 CN 








Jk = 5i - 4i 


91985.318 








CH3CN 








Jk = 52 - 4 2 


91979.998 








CH 3 CN 








Jk = 5 3 - 4 3 


91971.134 








CH 3 CN 








Jk = 60 - 5 


110383.504 




30£ 


32 d 


CH3CN 








Jk = 61 — 5i 


110381.376 








CH 3 CN 








Jk = 62 - 52 


110374.993 








CH3CN 








Jk = 63-53 


110364.358 








HNCO 








JK a ,K c = 4(1,4 — 3o,3 


87925.237 




< 1.8 a 


25 


HNCO 








JK a ,K c = 5o,5 — 4o,4 


109905.749 




12 


23 


13 CS 








J = 2 — 1 


92494.308 




36 


32 


C 34 S 








J = 2 - 1 


96412.950 




73 


61 


C 33 S 








J = 2 - 1 


97172.064 




21 


14 


cs 








J = 2 — 1 


97980.953 




290 


210 


SiO 








J = 2 — 1 


86846.960 




13 


41 


NS 








2 n 1/2 N=2-l, J=5/2-3/2 


115153.835 




8.5 


17 



F=7/2-5/2, c-state 



a 3cr-upper limit 

b blend with C 2 H N = 1 - 0, J = 1/2 - 1/2, F = Q - 1 at 87407.165 MHz 
c Integrated emission for all J = 5 — 4 transitions 
^ Integrated emission for all J = 6 — 5 transitions 
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Table 1 - continued 



Molecule 


Transition 


Frequency 
(MHz) 


Methanol Integrated Intensity 
Maser (Kkm/s) 
Class NGC6334I NGC6334I(N) 


H 2 CS 


JK a ,K c = 3i, 3 - 2i,2 


101477.810 


>24 


>20 


H 2 CS 


JK a ,K a = 3l,2 — 2i : i 


104617.040 


27 


34 


SO 


Nj = 2 2 - li 


86093.950 


22 


19 


SO 


Nj = 2 3 - 1 2 


99299.870 


53 


69 


so 


Nj = 3 2 - 2i 


109252.220 


18 


21 


HC 3 N 


J = 10 - 9 


90979.023 


38 


48 


HC 3 N 


J = 11 — 10 


100076.392 


27 


38 


HC 3 N 


J = 12 — 11 


109173.634 


33 


59 


OCS 


J = 7 - 6 


85139.103 


12 


17 


OCS 


J = 8 - 7 


97301.209 


19 


19 


OCS 


J = 9 - 8 


109463.063 


19 


25 


CH 3 OCHO e 


JK a ,K c = 100,10 - 90,9 


111171.636 


4.5 


< 0.9 a 



e Blend of A and E types 



region (eg. I Walsh et ai]|l99Sft . it is clear that NGC6334I 
appears to be more evolved than I(N): it contains a promi- 
nent UC Hn region, as well as a strong infrared source , 
whereas only weak rad io continuum dRodriguez et al"1l2007l ) 
and infrared emission (IWalsh et al.ll 1999 ) has been found in 
the vicitiy of I(N). We also find evidence for an evolutionary 
difference in our data. Both the CS/N 2 H+ and CS/HNCO 
ratios are much larger for NGC 633 4 1 than for I(N). Previous 
work on the CS/N2H + ratio llZinchenko. Caselli fc Pirogovl 
l2009h and CS/HNCO ratio i|Martm et alj|2008h have indi- 
cated that CS is likely to be enhanced in abundance to- 
wards photodissociation regions (PDRs), whilst both N2H + 
and HNCO are likely to be reduced in abundance towards 
PDRs. Thus, NGC 6334 1 shows PDR-like ratios, presum- 
ably due to the UC Hn region, whereas abundances of both 
N 2 H + and HNCO do not appear to be reduced, compared 
to CS in NGC 6334 1(N), where a PDR is yet to form. 

When comparing the non-masing lines of CH 3 OH, we 
find that the emission is stronger towards NGC 6334 1 in the 
100.638, 107.159 and 111.289 GHz transitions. CH 3 OH emis- 
sion is stronger towards NGC 6334 1(N) in the 95.915, 96.741 
(quadruplet) and 97.582 GHz transitions. The stronger tran- 
sitions in NGC 6334 1 all have upper energy levels be- 
tween 103 and 304 K, whilst all the stronger transitions 
in NGC 6334 1 (N) have upper energy levels no greater 
than 22 K. Thus, we can attribute the stronger emission 
in NGC 6334 1 to hot CH 3 OH and the stronger emission 
in NGC 6334 1(N) to cold CH 3 OH emission. Whilst pre- 
vious work has shown that both NGC 6334 1 and I(N) 
have hot component s, with temperatures in excess of 400 K 
(|Beuther et alj|2007h . it is clear that the large-scale CH 3 OH 
observed here has only had time to heat up in NGC 6334 1. 

We detect CH 3 OCHO only towards NGC 6334 1. This 
molecule is thought to be produced only after a young star or 
protost ar has heated up the su rrounding material to above 
100 K (|Garrod fc HerbstJ I2006T ). This provides further evi- 
dence of the more evolved state of NGC 6334 1 over that of 
I(N). 



4.3 Line ratio maps 

The ratio of integrated intensities for different spectral lines 
can potentially be used to identify physical and chemical 
characteristics of regions within the map. For example, a 
simple comparison of ratios for different lines of the same 
species might be used to map physical quantities such as 
density and temperature. However, most species detected in 
multiple spectral lines (ie. CN, CH 3 CN, HNCO, SO, HC 3 N 
and OCS) show only a narrow range of energy levels over 
the multiple spectral lines. Thus only small changes in line 
ratios are expected from different spectral lines of the same 
species. In addition to this, small differences in the relative 
position of maps made at different frequencies can lead to 
artifacts that can dominate line ratios. In the present study, 
these complicating factors make such line ratio maps too 
unreliable to interpret. The only exception to this is CH 3 OH, 
where clear differences in the intensities of lines are seen 
between NGC 6334 1 and I(N), as discussed in $4~T\ 

It is possible to compare line ratios between species, 
that have been observed simultaneously, as this eliminates 
any relative positional offsets, as well as any varying effects 
of weather. In Figure [4j we show images based on the ratio 
of line integrated intensities for different species. 

The upper-left image shows the C ls O/ 12 CO line 
ratio map. The greatest ratio is found surrounding 
NGC 6334 1(N). A high ratio between these two lines means 
that the 12 CO line is optically thick. We also see a secondary 
peak in the ratio to the south-west of NGC 6334 E, which is 
close to the peak 12 CO peak integrated intensity (shown in 
contours). At this position, it is likely that there is a high 
column density of gas. At NGC 6334 1, the line ratio is rela- 
tively low. The 12 CO gas is likely optically thick here since 
it is one of the centres of high column density gas in many 
other lines. The relatively low ratio, compared to that found 
in NGC 6334 1(N) indicates that the 12 CO emission is likely 
to be optically thick along most lines of sight cvered in these 
observations. 

The upper-right image shows the HNCO/C 18 line ra- 
tio map. As mentioned in §4.21 HNCO tends to be destroyed 
in PDRs. Thus, the line ratio map broadley shows two re- 
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Figure 4. Line ratio maps. Upper-left shows the C ls O (1— 0)/ 12 CO (1—0) ratio in greyscale and 12 CO (1-0) integrated intensity in 
contours. Upper-right shows the HNCO (5o,5— 4o,4)/C ls O (1—0) ratio in greyscale and C ls O (1-0) integrated intensity in contours. 
Lower-left shows the CH3CN (5-4)/N2H+ (1-0) ratio in greyscale and N2H+ integrated intensity in contours. Lower-right shows the 
13 CS (2-l)/N2H+ (1—0) ratio in greyscale and N2H+ (1-0) integrated intensity in contours. Dark areas show high line ratios and light 
areas show low line ratios. The same contour levels as for Figure[3]have been used. Each ratio map is masked, based on the weaker image 
(the numerator), to reduce contributions due to non-real emission. However, non-real emission is seen at the bottom of the line ratio 
map for CH3CN/N2H+ map, below the N2H+ contours. 



gions in the map: a region to the east of NGC 6334 1(N) 
which shows a high ratio (ie. not a PDR) and a region 
including NGC 6334 1 and E, which appears to be a good 
candidate PDR, although low signal in HNCO limits our in- 
terpretation beyond a narrow strip joining NGC 6334 1 and 
E. 

The lower-left image shows the CH3CN/N2H + line ratio 
map. As mentioned in H4.ll the CH3CN emission is strongly 
confined to the centres of active star formation: NGC 6334 1 
and I(N), whereas N2H + tends to trace ubiquitous and qui- 
escent gas. The line ratio map clearly shows low line ra- 
tio towards NGC 6334 1, in contrast to a high ratio seen at 
NGC 6334 1(N). This demonstrates that I(N) is a more qui- 
escent region, presumably with star formation at a younger 
stage. 

The lower-right image shows the 13 CS/N2H + line ratio 
map. Here, we see the highest line ratio at the position of 



the Hn region NGC 6334 E. This is probably the result of 
N2H + being destroyed by the Hn region. 



5 CONCLUSIONS 

We have mapped a 5' x 5' region, encompassing the high 
mass star formation sites NGC 6334 1 and I(N). We have 
covered most of the 3 mm spectral window with these obser- 
vations, using the Mopra radiotelescope. Our observations 
do not have sufficient spectral resolution to resolve spec- 
tral lines, but we compare morphologies of the various lines 
detected. We detect emission from 19 different moleules, 
ions and radicals, including multiple transitions of the same 
species and transitions from various isotopologues to yield a 
total of 52 spectral line detections. 

We find that CH3CN most closely follows the sites of 
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active star formation, although these sites feature promi- 
nently most other observed species. In contrast, emission 
from both CN and C2H appear to be widespread and are 
good tracers of gas that is not associated with active star 
formation. N2H + and HC3N morphologies closely resemble 
that of dust continuum emission, indicating that these two 
species are good tracers of column density, associated with 
both star forming sites and quiescent gas. 

Interpretation of the emission from common species and 
their isotopologues (ie. CO, HCO+, HCN and HNC) are 
difficult to interpret, which is likely due to some transi- 
tions being optically thick. All CO isotopologues show a 
peak of emission about 45" to the south-west of the Hn 
region NGC6334E. However, this peak is not prominent in 
any other tracer, including dust continuum emission. We are 
currently unable to provide an explanation for this. 

CH3OH emission is seen in many transitions, includ- 
ing known Class I and II masers. Due to the low spatial 
resolution of these observations, we cannot conclude that 
we have detected any new masers amongst these new tran- 
sitions. Thermal CH3OH transitions show a marked dif- 
ference in their ocurrence: transitions with upper level ki- 
netic temperatures above 100 K are found to be stronger 
towards NGC 6334 1 and transitions with upper level ki- 
netic temperatures below 22 K are found to be stronger to- 
wards NGC 6334 1(N). This supports the interpretation that 
NGC 6334 1 is a more evolved site for star formation than 
I(N). 
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